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Abstract 

We calculate e^e~ QQ{k)-\- anything in a certain momentum, k, region 
for heavy quark-antiquark (QQ) production. In our chosen region we find 
that the number of heavy quark pairs produced is determined by BFKL 
dynamics and the energy dependence of the number of pairs is given by 
ap — 1 the hard pomeron intercept. 

1. Introduction 

The topic of this paper, BFKL dynamics [1, 2] in QCD jet evolution 
[3, 4] is perhaps surprising since it is generally thought that BFKL dynamics 
has uniquely to do with high-energy scattering and not at all to do with 
jet physics. Indeed, to a large extent jet evolution is an understood and 
well-tested subject with the dominant dynamics being given by double log- 
arithmic perturbative terms coming from emission off the primary energetic 
partons. Their resummation leads to a Sudakov type of distribution. Until 
very recently there has been no hint of BFKL dynamics in jet physics. How- 
ever, about two years ago Dasgupta and Salam [5] observed that in order to 
calculate certain "non-global" observables such as the probability, S, that a 
jet of energy E decay while emitting less than a certain amount of energy, 
Eo^t, into a specified angular region away from the jet, it is necessary to 

-'^This research is supported in part by the US Department of Energy 



1 



evaluate single logarithmic terms in the jet evolution. For S double logarith- 
mic terms, soft and coUinear, cancel leaving only the logarithms from soft 
non-coUinear emissions originating from successive branching of soft gluons 
at comparable angles. 

Stimulated by the results of Dasgupta and Salam, Banfi, Marchesini and 
Smye (BMS) [6] found an equation for all the single logarithmic terms con- 
tributing to S derived from the multi-soft gluon distributions given in Rcf. [3]. 
The surprise is that the dominant part of the BMS equation is identical to 
the Kovchegov equation [7]. This is surprising for several reasons: (i) The 
Kovchegov equation is an approximate equation for dealing with small-x 
evolution, or high-energy scattering, when unitarity limits are being reached. 
The BMS equation is not an approximate equation; it is an equation which 
exactly deals with single logarithmic terms. However, it appears to have 
nothing to do with unitarity limits, rather, the nonlinear aspects of the equa- 
tion are related to the nature of the non-global observable S. It is due to the 
non linearity of successive branchings of soft gluons (or jets) at comparable 
angles, (ii) The linear version of the Kovchegov equation is the BFKL equa- 
tion, and when scattering is far from unitarity limits the Kovchegov equation 
becomes the BFKL equation. There appears to be no limit in which the BMS 
equation reduces to the BFKL equation. The linear limit of the BMS equa- 
tion is just low orders of perturbation theory. 

Is it just accidental that the BMS equation is the same as an equation 
which naturally appears in high-energy scattering, or is it rather an indication 
that BFKL dynamics is also present in jet decays, waiting to be found? We 
believe the latter is the case, and we have found an observable in jet decays 
which, in certain kinematic regions, is dominated by BFKL dynamics. The 
observable is easiest to describe in the context of e^e~ annihilation, although 
we could equally well define similar observables in deep inelastic scattering 
or in hadron-hadron collision events with large Et- In the center of mass of 
the collision let Eq be the energy of the e+e" pair. The observable is the 
number of heavy quark-antiquark pairs produced having pair mass M. on the 
order of 2M, with M the heavy quark mass, and with /cq, the energy of the 
pair, in the regime ko/Ad on the order of one. We suppose that Eo/ko » 1. 

Of course most heavy quarks are not produced in the kinematic region 
where we are looking. The majority of heavy quarks will be collinear with 
one of the two primary jets and their production will be dominated by double 
logarithmic terms. Our choice of kinematics is motivated by the fact that 
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double logarithmic terms cancel out in the region we have taken. In subse- 
quent sections, where the calculation is performed, we choose, for technical 
reasons, a boosted frame where the small angle approximation can be used. 

In section 2, we calculate the lowest order contribution of a heavy quark- 
antiquark emission from a QCD dipole. 

In Section 3, we calculate the evolution of the two-jet system, originally 
produced in the e^e~ collision, in terms of QCD dipoles. This section makes 
use of the large Nc approximation, however, it may well be that similar results 
can be obtained without the large Nc approximation. 

In section 4 we compare the two-jet evolution equation with the dipole 
version of QCD evolution in high-energy scattering and find them to be 
formally equivalent, but in terms of different variables. In jet evolution it is 
the angles of gluons which are frozen during subsequent emissions and hence 
become the natural variables, while in high-energy scattering it is transverse 
coordinates which are frozen and so are the natural variables. 

In section 5 we give the rate of heavy quark- ant iquark pair production 
in our chosen kinematic region, and we find that ap — 1, the hard pomeron 
intercept, determines the jet energy dependence of the production. 

2. The observable and the Born approximation 

The observable we consider is heavy quark-antiquark production in jet 
decay. The calculation will be done in the large Nc (planar) approximation 
where one can view jet evolution in terms of a sequence of softer and softer 
gluon emissions which in turn can be viewed as the production of more and 
more QCD dipoles. It may be that much of what we derive can also be 
obtained beyond the large Nc limit, but that is beyond the scope of the 
present paper. The simplicitiy of the large Nc limit is that one may view the 
gluon as being made of a quark part and an antiquark part. 

In this section we focus on the heavy quark-antiquark pair (QQ) emission 
from a dipole consisting of the quark part of a gluon of momentum Pa and 
the antiquark part of a gluon of momentum pij. The gluons Pa and pb have 
come from earlier parts of the decay of a jet, or pair of jets, and it is assumed 
that IPaMPbl ^ 1^1 = with the four-momentum and ^^/k^ki^ — M. the 
mass of the QQ pair. We denote by pi and p2 the momenta of the Q and 
Q, respectively, with k^ = {pi +P2)n- It is important to emphasize that k is 
the softest momentum in our problem, softer than all other gluons emitted 
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in the secondary branchings, so below we are able to use the soft momentum 
approximation. 



We now turn to the calculation of the lowest order contribution to pair 

production from the dipole {pa,Pb)- We denote this contribution by ^j^^l'^y , 
where y — Iuk/M. One of the lowest order graphs is shown in Fig. 1 where 
the left-hand part of the graph denotes the amplitude and the right-hand 
part of the graph the complex conjugate part of the amplitude. The factors 
accompanying the coupling of the gluon, labelled by k, to the "quark" a or 
to the "antiquark" b are given by 

{2pa + k)a Paa {2pb+k)a Pba v 

[Pa + k)^ Pa -k {pb + ky Pb ■ k 

with color factors to be put in later. In writing (1) we assume a leading 
soft approximation where |/c|/|pa| and ^ 1, but we do not assume 

coUinear emission. Indeed, we shall in the end fix AI/k = so as to 

guarantee that coUinear emission is unimportant. 
One can write 



dN^^^ r d^k 



dM'^dy 



Paa Pba \ / Pad Pb(3 



Pa-k Pb - kj \pa - k Pb - k 



2Im7r°^^(A;) 



(2) 



where 

2Im7r„/3(A;) = 2 {k^gap - kakp) Im7r(/c^) 

= ^/ ||^<^'(^-Pi-P2)tr{(7-Pi+M)7.(7-P2-M)7,} 
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with M the heavy quark mass. From (3) it is straightforward to get [8] 



Im7r(A^^) = 



(4) 



Because of current conservation the kakp term in Im7rQ^(/c) does not con- 
tribute to (2) leaving 



dN' 



(0) 



ah 



dM^dy 



ATrasCp 



d^k 



(27r)47W4 



6 {k,k^ - M') S{y-\n 



(5) 



{Pa ■ k){k-pb) 
It is now staightforward to get 



dL dQk 



I 



(6) 



where 6ak is the angle between the vectors Pa and k, with similar definitions 
for 9kb and 9ab, and where we take M./n <^ 1, and we suppose 6ah <^ 1 so 
that the small angle approximation cos^ab ~ 1 — ^^^^ along with similar ap- 
proximations for cos Oak and cos 9kh can be used. Finally, it is straightforward 
to evaluate 



/ 



el,dVtk ^ 27rA 

This gives 



In 



A 



ab 



2M 



(7) 



dN, 



(0) 



ab 



dM^dy 37r2 



alCp VM^-AM^ ( 2M2 



A 



) ^/T^:p 



In 



Tl+F+A 
a/T+A2-A 



(8) 



We imagine choosing TW^— 4M^ on the order of so that the counting of 
heavy quarks and heavy mesons should be the same and physically observ- 
able. We also choose A on the order of one so that there are no collinear 
singularities in the emission of the heavy quark pair from the dipoles. 
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3. The evolution equation 

The evolution (branching) [3, 9] of QCD jets has been widely studied and 
furnishes some of the best tests of perturbative QCD. Double logarithmic 
terms, having both soft and coUinear singularities, dominate the behavior of 
most global observables, and the resummation of these perturbation series is 
now well understood. 

The evolution we are going to review in this section is single logarithmic. 
Soft, but not coUinear singularities, will be summed. Such terms govern 
the behavior of certain non-global observables such as -Eout, the total energy 
emitted into a region Cout away from all the hard jets. The same single 
logarithmic terms will dominate heavy quark-antiquark production in certain 
regions of phase space as we shall see below. However, for the moment we 
are going to describe more formally the evolution that will be dominant for 
our heavy quark-antiquark production. We shall see that the evolution is 
determined by the formalism developed in by Banfi, Marchesini and Smye 
(BMS) [6] in their discussion of the distribution function for Eout mentioned 
above. We shall also see, perhaps surprisingly, that the evolution is formally 
exactly equivalent to BFKL evolution as described in the dipole formalism. 

We begin by considering the process e'^ + e~ 1*{q) ~^ quark {pa) + 
antiquark (p^) . We chose an unusual frame where q is large and along the 
z-axis. Without loss of generality we may suppose that \pa\ — \pb\ — E for 
our zero mass quarks and that 

- £; E 

is small. Q, the largest physical scale, is assumed to be much greater than 
the heavy quark mass M introduced in the previous section. 

Of course, the quark-antiquark pair {pa,Pb) will further evolve by gluon 
emission. The gluon emission probabilities are known analytically [3] in 
the large Nc limit and in the leading soft approximation for longitudinal 
momenta. No assumption of a coUinear approximation is necessary. In the 
large Nc approximation the states having n gluons in addition to the original 
quark-antiquark pair may be viewed as a state of n -|- 1 dipoles. For example, 
the state having a quark-antiquark pair and a single gluon consists of two 
dipoles. The first dipole is the original quark and the antiquark part of 
the gluon while the second dipole consists of the quark part of the gluon 
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and the original antiquark. Since our object is to calculate heavy quark- 
antiquark production, and since the heavy quark pair can only be emitted 
from a single dipole, the same dipole in the amplitude and in the complex 
conjugate amplitude, what we need to extract from the QCD evolution is 
only the inclusive dipole distribution caused by the gluon emissions. 

Therefore, we need to compute n{9ab, d, Y) the mimbcr density of dipoles 
of opening angle 6, starting from the original quark- antiquark pair which has 
opening angle 9ab, and in a rapidity interval Y — In E/k. (Here we take 
the rapidity of a quark or gluon of energy a; to be y = lna;/A.) The dipole 
having opening angle 6 consists of the quark part of a gluon, say gi, and 
the antiquark part of a gluon, g2. (The original quark or antiquark may also 
make up either the quark or antiquark part of the dipole 6, but this is unlikely 
to be the case at large Y. Our formahsm allows this possibihty.) Let cv be 
the smallest of the energies of the gluons gi and g2- Then in n{9ab, 0, Y) we 
require that hiE/uj <^ Y. The result for the distribution of heavy quark pairs 
is 

J A, 



dMW -J ^ ^^''''^^ dMMy ' 

where y — Ihk/ M. and the details of our normalization for n will be explained 
later. 

In order to write an equation for n we need only an explicit expression 
for the emission of a single gluon, of momentum A:, from the original quark- 
antiquark pair. By a simple calculation (see for example, Eqs. 3.2-3.4 of 
BMS) , one finds for the number of gluons Ng 

_ , dVtk Pa-Pb - ^sNc .^„s 

dNg^a^udu^j—^^^^, o^s--^. (10) 

(In covariant gauge only the two graphs shown in Fig. 2 contribute.) Going 
to angular variables 

(11) 

in the small angle approxiation, and where y = In a;/ A. 

Eq. (11) can be turned into an equation for the inclusive dipole density 
simply by observing that the measured dipole has as a parent dipole either 



7 



a + antiquark part of k or quark part oi k + b. Thus, as illustrated in Fig. 3, 



dn{6ab, 9, Y) _ 
dY ~ 2^ 



[niOak. 0, Y) +n{9k,, 9, Y) -n{9ab, 9, Y)] , (12) 



where the third term on the righthand side of (12) is the virtual contribution 
of the gluon k. Eq. (12) is our basic equation, and we shall come back to 
it when calculating an explicit asymptotic result for heavy quark- antiquark 
production. 




FIG. 2a 





Fig. 3 



In Fig. 3 the double line, k, represents the gluon, k, emitted coherently from 
a and b. The final term on the righthand side of Fig. 3 has the k-line being 
emitted and then reabsorbed so that it does not appear in the final state. 

4. Relationship to BFKL dynamics 
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Eq. (12) bears a remarkable resemblance to the BFKL equation as given 
in the dipole formulation [10, 11, 12]. The BFKL equation reads 



dn{xoi,x,Y) as fxQid X2 

'^02^21 



dY 2n. 



— — [n{xo2,x,Y)+n{x2i,x,Y)-n{xoi,x,Y)], (13) 

^02'*^21 



for the number density of dipoles of size x to be found in a parent dipole of 
size Xqi over a rapidity range Y, where the parent dipole is constructed from 
a quark at Xq and an antiquark at Xj^ with xqi = |xq — X]^|. The measured 
dipole has a transverse coordinate separation x between its quark and anti- 
quark parts. In (13) we anticipate that the inclusive dipole distribution is 
independent of the orientations of the dipoles when Y is large. (In (12) we 
have similarly anticipated that n only depends on the magnitudes 9ab and 
6, but not on the orientation between the angles of the dipoles.) From (13) 
it is straightforward to evaluate high-energy dipole-dipole scattering whose 
rapidity (energy) dependence is the same as n in (13). 

Except for the measure of integration dflk versus d'^X2 equations (12) and 
(13) are identical. However, in the small angle limit, 

dQk - d(l)kOkdek (14) 

can be viewed as an integration measure over a fiat plane with 9^ and (pk the 
radial and polar coordinates respectively. Thus (12) and (13) are formally 
identical equations so long as we allow the ^jt-integration to cover the region 

< ^fe < oo (15) 

in (12) while still using (14). Of course when 6ah is small the important 
region of integration in (12) is 9k < 9ab so that (15) is necessary only to make 
(12) and (13) exactly the same. 

Although (12) and (13) are formally indentical the physical variables 
which appear in these two equations are very different. In jet evolution angles 
seem to be the preferred variables while in high-energy scattering transverse 
coordinates seem to be more natural. It is not hard to see why this is the 
case. In both cases the basic amplitude is the emission of a softer gluon off 
a higher momentum quark or gluon as illustrated in Fig. 4 for a momentun 
labelling of these gluons. 
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p + k 




Fig. 4 



In the high-energy scattering case the graph in Fig. 4 represents a gluon (p) 
+ gluon (k) part of a hght-cone wavefuncton. The soft emission hmit means 
k+/{p + k)+ = z <^ 1. The transverse momenta k and p are typically of the 
same size, and in the double logarithmic region + "C |^|. The light-cone 
quantization energy denominator is 



D~^ = \p_ + k-- (p+k). 



p 



2p+ 2k+ 2(p+k) 



+ . 



-, -1 


\- 1 


-1 




_2k+_ 





(16) 



and is dominated by the softest gluon, k. Thus, the time over which the 
gluon, k, is being emitted is 

r.-^. (17) 

The transverse velocity of the gluon p is = p/p+ while that of the gluon k 
is Vk = k/k+- Over the time during which k is being emitted the changes in 
the transverse coordinate positions of the p and /c-lines are 



and 



I Ax I ~ \v\ Tfe OC — 1^1 ^ 



\AxJ ~ 1^1,1 Tk oc \k\ 



(18) 



(19) 

respectively. Eq. (19) reflects the uncertainity principle while (18) says that 
the harder gluon has a very small change in its transverse coordinate during 
the emission of softer gluons. It is thus a good approximation to "freeze" 
the transverse coordinate of a high momentum quark or gluon in a light- 
cone wavefunction during the time of formation of the softer parts of the 
wavef unction. 

On the other hand in jet evolution the energy denominator is better 
written as 

z(l-z) 

- I 7)H- k: I 1 I 



D = 



{p+k). 



(20) 
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with 

"~{p+kh' "~{p+kh ^ ^ 

and 

P V2p V2k 

Pz P+ 

In jet evolution 9_p and 9_j, are typically of the same size, with \ 9_p — 9_|^^ <^ \9_p\ 
in the double logarithmic limit so that typically |^| oc z \p\ for soft gluon 
emission. This means that 



^p = - = -r^' (22) 



{p + k)+ ~ p+ 



(23) 



for soft gluon emission. Eq. (23) indicates that the angle of a harder gluon 
does not change due to softer gluon emissions and so angular variables nat- 
urally appear in the jet evolution equation (12). 

One final issue in comparing high-energy evolution with jet evolution is 
the question of the running of the QCD coupling. Equations (12) and (13) 
have been written for fixed coupling evolution. Running coupling effects come 
in very differently in (12) and (13). In (12), following BMS, we introduce a 
variable A where 

du , 2^2N d Ed Id 

and (12) becomes (replacing n{9ab, d, Y) n{9ab, d, A)) 



dn{0„,.9.A) fd^k el, 



dA J 2n ei,ei 



n{9ak, 9, A)+ni9kb, 9, A)-n(^„,, ^, A) ] . (25) 



When 9 and 9ab are of similar magnitude (24) and (25) should be adequate 
to represent running couphng effects. We suppose, of course, that E > Eq 
and that EQ9ah/ ^qcd ^ 1- (Note that in jet evolution the diffusion of angles 
away from 9 and 9ab is not a serious problem in determining the argument of 
the running coupling since the cu-dependence in (24) dominates the variation 
of a during the evolution). Thus, in jet evolution running coupling effects 
are included, when 9ab and 9 are not too different, simply by the replacement 

cXsY^A 
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in going from the fixed coupling to the running coupling case. 

Running coupling effects in high-energy BFKL applications are reason- 
ably well understood [13, 14, 15]. For our purposes it is sufficient to ob- 
serve that for Xqi and x of comparblc size one may take ~ aJ^(xQi) ~ 
— 61n(a;Q^AQ(^p) so long as a^Y^ <^ 1. When a^Y^' ^ 1 high-energy evolu- 
tion and jet cvohition will look quite different. We note that unitarity effects 
become important in the regime a^Y^ <^ 1. 

5. The leading soft approximation result 

It now only remains to combine our lowest order result for heavy quark- 
antiquark production from a single dipole with the solution of (12) for large 
values of Y. For large Y the solution of (12) is given by the standard BFKL 
formula [12] 

n{e^,, 0,y)-l,j J^n. e^^-M^ l^-f^ (26) 

where 

X{v) = m - h^ih - - ¥i\ + (27) 
and flu determinines the normalization of n at F = 0. If we take = 4 then 



n{eab, 0,0) (\n Y-)=l ^(^-M , (28) 

in which case, in a notation and normalization used in (9), 

J ^n(^„,,^„,,,,F = 0) = l. (29) 
It is straightforward to find 

with, as usual, 

ap - 1 = ln2. (31) 

In case of running coupling evolution, Eq. (25), one simply replaces (isY by 
A in (30). However, in the running coupling case it is difficult to get large 
values of A in realistic circumstances. 
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Now using (8) and (30) in (9) one arrives at the number of heavy quark 
pairs in the original dipole to be 



dNab 



) 



dM^dY 



24: \M 



|«,7VeC(3)y 



where Y = In E/k, and we remind the reader that Ai is the mass of the 
heavy quark- antiquark pair and k = \k\ the pair momentum and we suppose 
that nOab/ M. is not too different from 1 in order to stay within the (angular) 
diffusion radius for BFKL evolution. 
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